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Abstract. The location and asymptotic behaviour for large n of the zeros of exceptional Jacobi 
and Laguerre polynomials are discussed. The zeros of exceptional polynomials fall into two 
classes: the regular zeros, which lie in the interval of orthogonality and the exceptional zeros, 
which lie outside that interval. We show that the regular zeros have two interlacing properties: 
one is the natural interlacing between consecutive polynomials as a consequence of their Sturm- 
Liouville character, while the other one shows interlacing between the zeros of exceptional and 
classical polynomials. A generalization of the classical Heine-Mehler formula is provided for the 
exceptional polynomials, which allows to derive the asymptotic behaviour of their regular zeros. 
We also describe the location and the asymptotic behaviour of the exceptional zeros, which 
converge for large n to fixed values. 



1. Introduction 

Let n be a probability measure supported on an infinite subset E of the real line. We assume that 
J E \x n \dfi < oo for every nonnegative number n. A sequence of monic polynomials {P n (x)}n>a is 
said to be orthogonal with respect to fi when deg P n = n and J E P n (x)P m (x)d^i(x) = for n ^ m. 
It is very well known that these polynomials satisfy a three term recurrence relation that yields for 
the orthonormalized polynomials a symmetric tridiagonal (Jacobi) matrix such that the eigenval- 
ues of the n leading principal submatrix are the zeros of the polynomial P n . As a straightforward 
consequence of this fact the zeros of P n are real, simple and interlace with the zeros of P n -i- On 
the other hand, they are located in the interior of the convex hull of E. 



The theory of orthogonal polynomials is strongly related to Sturm-Liouville problems. In par- 
ticular, the so called classical orthogonal polynomials (Hcrmite, Laguerre, and Jacobi) appear as 
eigenfunctions of second order linear differential operators with polynomial coefficients. Indeed, the 
corresponding measure of orthogonality is absolutely continuous and their derivative with respect 
to the Lebesgue measure (weight function) is the density function of the normal, gamma and beta 
distributions, respectively. Notice that this fact was pointed out by E. Routh in 1884 41 as well 
as by S. Bochner in 1929 [2], but the orthogonality does not play therein any role. On the other 
hand, they are hypergeometric functions and, as a consequence, many analytic properties can be 
deduced from this fact. Moreover, certain properties of their zeros can be easily deduced using 
the classical Sturm theorems. Finally, a nice electrostatic interpretation of their zeros is deduced 
from the second order linear differential equation as an equilibrium problem for the logarithmic 
interaction of positive unit charges under an external field. 



Exceptional orthogonal polynomials constitute a recent new approach to spectral problems for 
second order linear differential operators with polynomial eigenfunctions. Previously, a construc- 
tive theory of orthogonal polynomials related to the classical ones has been done in two directions. 
The first one is related to the spectral theory of higher order linear differential operators with poly- 
nomial coefficients. For fourth order differential operators the classification of their eigenfunctions, 
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which are sequences of orthogonal polynomials with respect to a nontrivial probability measure 
supported on an infinite subset of the real line, was done by H. L. Krall and A. M. Krall [27 29 
and essentially yields the classical ones and perturbations of some particular Laguerre weights 
c - * + MS(x), Jacobi weights (1 — x) a + M5(x) and Legendre weight 1 + MS(x — 1) + MS(x + 1), 
M > 0. For higher order, some examples are known but a general theory and classification con- 
stitutes an open problem. The second one appears when some perturbations of the measure are 
considered. In particular, three cases are considered in the literature in the framework of the so 
called spectral linear transformations 45 . The Christoffel transformation (the multiplication of 



the measure by a positive polynomial in the support of the measure), the Uvarov transformation 
(the addition of mass points off the support of the measure) and Geronimus transformation (the 
multiplication by the inverse of a positive polynomial). They can be analyzed in terms of the 
discrete Darboux transformation of the corresponding Jacobi matrices using the LU and UL fac- 
torizations and commuting them (3l. 



Exceptional orthogonal polynomials depart from the classical families in that the sequence of 
exceptional polynomials is not required to contain a polynomial of every degree, and as a conse- 
quence new differential operators exist, with rational rather than polynomial coefficients. Despite 
this fact, the sequence of exceptional polynomial eigenfunctions is still dense in the corresponding 
weighted L 2 space and constitutes an orthogonal polynomial system. The measure of orthogonality 
for the exceptional families is a classical measure divided by the square of a polynomial with zeros 
outside the support of the measure. 

The first explicit examples of families of exceptional orthogonal polynomials are the Xi-Jacobi 



and Xi-Laguerre polynomials, which are of codimension one, and were first introduced in 15 



16 



In these papers, a characterization theorem was proved for these orthogonal polynomial 
families, realizing them as the unique complete codimension one families defined by a Sturm- 
Liouville problem. One of the key steps in the proof was the determination of normal forms for the 
flags of univariate polynomials of codimension one in the space of all such polynomials, and the 



determination of the second-order linear differential operators which preserve these flags 13 19 



Shortly after, Quesne 33 34 observed the presence of a relationship between exceptional or- 
thogonal polynomials and the Darboux transformation^] This enabled her to obtain examples 
of potentials corresponding to orthogonal polynomial families of codimension two, as well as ex- 
plicit families of X2 polynomials. Higher-codimensional families were first obtained by Odake and 
Sasaki 36 . The same authors further showed the existence of two families of X m -Laguerre and 
X m - Jacobi polynomials 37 , the existence of which was explained in 17 for X m -Laguerre polyno- 
for X m -Jacobi polynomials, through the application of the isospectral algebraic 

These exceptional orthogonal polynomials 



19 



Darboux transformation first introduced in 11 12 



have been applied in a number of interesting physical contexts, such as Dirac operators minimally 
coupled to external fields, 24 , entropy measures in quantum information theory, |9j, rational 



21 22 or discrete quantum mechanics, 40 



extensions of Morse and Kepler- Coulomb problems, 

The aim of our contribution is to explore analytic properties of these exceptional polynomials. 
In particular we will focus our attention in the distribution of their zeros in terms of the support of 
the orthogonality measure as well as their limit behavior. On the other hand, we will analyze some 
asymptotic properties as the outer relative asymptotics in terms of the corresponding classical 
orthogonal polynomials and the Mehler-Heine type formulas. Some properties of the zeros have 



also been analyzed numerically in 25 . 



By Darboux transformation, wc do not mean here the factorization of Jacobi matrices into upper triangular 
and lower triangular matrices mentioned above, but the factorization of the second order linear differential operator 
into two first order linear differential operators 11 12 
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2. Exceptional orthogonal polynomials 

Let W(z) be a positive weight function with finite moments. Usually, orthogonal polynomials are 
defined by applying Gram-Schmidt orthogonalization to the standard flag 1, z, z 2 , . . . relative to an 
L 2 inner product associated with the weight W . Moreover, if the resulting orthogonal polynomials 
are eigenfunctions of a Sturm-Liouville problem, we speak of classical orthogonal polynomials. By 
Bochner's theorem, the range of such polynomials is limited to the classical families of Hermite, 
Laguerre, and Jacobi (for positive weights) and Bessel (for signed weights). 

In order to go beyond the classical families, we consider orthogonal polynomials spanning a non- 
standard polynomial flag, say with a basis p m (z) , p m+ i(z) , . . . , where degpj — j. Once we drop 
the assumption that the OP sequence contains a polynomial of every degree, we obtain new classes 
of orthogonal polynomials defined by Sturm-Liouville problems, which are commonly referred to 
as exceptional orthogonal polynomials (XOPs)[^j 

In the last two years it has become clear that the Darboux transformation, appropriately gener- 
alized to the polynomial context, plays an essential part in the deliniation of XOPs. To wit, let V n 
denote the vector space of polynomials of degree < n, and consider a codimension m polynomial 
flag 

U = {U k }f =1 , C/fe-i c U k c V m+k -x . 

Furthermore, let 

(1) A[y] = b(z)(y' - w(z)y), B[y] =b(z)(y' - w(z)y) 
be first order linear differential operators with rational coefficients such that 

(2) Bp> fc -i] = U k , A[U k ]=V k - 1 , fc = I,2,... 

i.e. B maps the standard flag into the codimension m flag while A maps the codimension m 
into the standard flag. Note that equation ([2| implies that ker A = kcr B = and the Darboux 
transformation is isospectral. Next, consider the second order differential operators 

(3) T = AB, T = BA. 

By construction, T leaves invariant the standard polynomial flag, while T leaves invariant the 
codimension-m flag U. Hence, by Bochner's theorem, 

(4) T[y] = p{z)y" + q(z)y' + r(z)y, f[y] = p(z)y" + q(z)y' + r(z)y 

where p,q,r are polynomials with degp < 2,degg < l,degr = but where q,f are ingeneral 
rational functions. It is then assured that T and T have polynomial eigenfunctions [j . Let 
Vj-, j ^ 0> denote the polynomial eigenfunctions of T[y] and yj, j > m, denote the polynomial 
eigenfunctions of T[y]. Again, by construction we have the following intertwining relations 

(5) TA = AT, BT = TB, 
which mean that 

(6) B[yj] = PjVj+m, A[y j+m ] = a 3 y , j = 0, 1, 2, . . . , 

where Oj, fij are constants, i.e. operator A maps eigenfunctions of T into eigenfunctions of T while 
B does the opposite transformation. 
Furthermore, let 

(7) W(z) - 1 exp /\ W = 1 exp 

P J P P J P 



2 Note that the requirement that the degree sequence starts at m and contains every integer j > m is not essential 
either, although all the families treated in this paper belong to this class. There exist also XOPs where the degree 
sequence has gaps. They are related to state-adding Darboux transformations (as opposed to isospectral) and 
contain for instance the X -Hermite families, beside many others. 

■^Since no domains have been specified for A and B, the term eigenj unction is not meant in the strict spectral 
theoretic sense here, but rather as polynomial solutions to the eigenvalue equation. 
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be the solutions of the Pearson's equations 

(8) ( P wy = q w, ( P wy = q w. 

This means that T[y] is formally self-adjoint relative to W while T[y] is formally self-adjoint relative 
to W . Consequently, the eigenpolynomials yj arc formally orthogonal with respect to the weight 
W(z), while jjj are formally orthogonal relative to W(z). One can also show that 

J A[f]gWdz + J B[g]fWdz= boundary term; 

which means that operators A and — B are formally adjoint. By a careful choice of the flags, and 
by imposing appropriate boundary conditions one can construct examples where the above formal 
relations hold in the L 2 setting (i.e. boundary conditions such that the boundary terms vanish) 
and thereby obtain novel classes of exceptional orthogonal polynomials. 

In the present note we study asymptotic behaviour of XOPs of Laguerre and Jacobi types. 
As we show, in the interval of orthogonality the exceptional polynomials satisfy a variant of the 
classical Heine-Mehler formula. Outside the interval of orthogonality, the convergence picture is 
less clear. However, one can show that codimension m exceptional orthogonal polynomials possess 
m extra zeros outside the interval of orthogonality, which we shall denote as exceptional zeros. 
These exceptional zeros have well-defined convergence behaviour, and they converge to the zeros 
of some fixed classical orthogonal polynomial. 

3. Type I Exceptional Laguerre polynomials 

Let us illustrate the above discussion with the particular example of the so-called type I ex- 
ceptional Laguerre polynomials. Let (z) denote the classical Laguerre polynomial of degree n 
and 

C a [y] =zy" + (a + l-z)y' 
the classical Laguerre operator. Thus, y — is the unique polynomial solution of the equation 

£a[y] = — n Ui with the normalization y^ n \0) = (— 1)". 
An equivalent boundary condition is 



(9) #>(0)=( n ; a ). 
For a fixed non-negative integer m > 0, let us now define 

(10) Z a , m {z) = L%X-z), 

(11) A a, m [y] = t a , m y' - £ a +i, m y, 

(12) BUy] = V + (l + a) y > 

(13) Ll% +j = -A l a _ l m [L^] , j = 0, 1, 2, . . . 

(14) ^mM = £<x[y] +my- 2(log^ Q _i, m ) / {zy' + ay) . 



The following factorization relations follow from standard Laguerre identities: 

(15) C a =B I aim A I aim + a + m+l, 

(16) C l a>m = A l a _ lm B l a _ lm + a + m. 
Consider now the following codimension m polynomial subspace 

(17) U I atj = {feP j+m - 1 :€ a - 1 , m \{zf' + af)}, J -1,2,..., 
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where f\g means polynomial f(z) divides polynomial g(z). At the level of flags, the above factor- 
izations correspond to the following linear isomorphisms: 



(18) 



B 



The polynomials I L 



i 

a— l,m 

}oo 
7 = 



At 



V 



3-1 



I (a 



are known in the literature as the type I exceptional codimension 



m Laguerre polynomials (for short, type I X m -Laguerre) [17|36 . By construction, the -X" m -Laguerre 
polynomials have the following properties: 

• they span the flag J7* ;1 C U^ 2 C • • • 

• they satisfy the following second order linear differential equation: 



(19) 



(20) 



r 1 



L 



1(a) 



• they are orthogonal with respect to the weight 

_j ( ^ z a e~ z z a e 

£,a-l,m(z) 2 r (a-l) 



3>0, 



Wl m {z) := 



they are dense in the Hilbert space L 2 ([0, oo), W a m ] 



E [0,oo), 



Note that for a > the polynomial £a-i,m(-2) — £m ( — ^) m the denominator of the weight 
has its zeros on the negative real axis, and hence W\ m dz is a positive definite measure in R + , 

with well defined moments of all orders. As a result, for a > the set {ll°l}^ =m constitutes 
an orthogonal polynomial basis of L 2 ([0, oo), W a m ) . We observe that the last property of the 
above list does not follow by the algebraic construction and needs to be established by a separate 



argument. The interested reader is referred to 17 for a direct proof of the completeness of the 
X m -Laguerre families. 

Proposition 3.1. The type I X m -Laguerre polynomials l}m)i can be expressed in terms of classical 
associated Laguerre polynomials with the same parameter a as follows 



(21) 



-1(a) 
J m,m-\-j 



-iL 



(a) 
3-1' 



j>0. 



The above representation will be specially useful to discuss the asymptotic properties of the 
zeros of type I A m -Laguerre polynomials. It is clear from it that L l ^l l+ - has degree m + j. 
This representation is reminiscent of the expansions obtained by rational modifications of classical 
weights in the framework of spectral linear transformations (see [45])- However, it should be 
stressed that they are essentially different because in the case of exceptional polynomials, although 
there is a rational modification of the weight, we are not dealing with the standard flag. 

Proof. Using elementary identities, we have 

n— t fr( a ) _ r( Q ' — r( a_1 )) _i_ rW (c _e _c ) 

u SC*,mv ;j j— 1 3 ' — lvsa.m ?a,m-l ^a-l,mj 



(a-1) 



Therefore, we can re-express the type I X m - Laguerre polynomials as 



(a) 



-(a) 



-1(a) 
J m,m-\-j 



t r(a-l) 

t r(«-l) 
SQ,m 

-(a) 



-(«-!)' 



-(«) 



Ca— 1,771-^7* — 1 

(a) 



Sa,m^j" ^a.m-l-^j-i 



□ 



We are now ready to prove an interlacing result for the zeros of type I A m -Laguerre polynomials, 
but before let us recall the following classical identity 

(22) 4 a >(0) = (a + l)„/n! 
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where 

, , [x(x + 1) • • • (a; + n- 1) if n > 0, 

W» = i , n 
II if n = 0, 

is the usual Pochhammer symbol. Using the above representation we obtain an analogous expres- 
sion for the type I A" m -Laguerre polynomials: 

Proposition 3.2. for a > £/ie type / exceptional Laguerre polynomial L TO " m+:) (z) has j simple 
zeros in z £ (0, oo) and m simple zeros in z € (— oo, 0). The positive zeros o/L m Q m+J (z) are located 
between consecutive zeros of L " and with the smallest positive zero of £ m Q m +j( 2 located 

to the left of the smallest zero of L^ a \ The negative zeros o/i m Q m+J (z) are located between the 
consecutive zeros of £ a ,m-i an d £a,m- 

Proof. Let < Cn.i < Cn.2 < ■ • • < Cn.n denote the zeros of L [ n\z) listed in increasing order. 
According to the interlacing property of the zeros of classical orthogonal polynomials, we have 

££i < (7-1,1 < Cj"a < Cj-1,2 < • ■ ■ < Cj-ij-i < Qj- 

Recall that Ln(z) > for z < 0. This implies that £ Q m (^),^ a m-i( z ) > for z > 0. Hence, by 
(21). 

orrn J 



sgn^l + ,(C, c ! i ) = (-l) i , ' 1 ./: 



sgnL^ + ,(C a _ M ) = (-!)*, / 1 j 1. 



m,m-\-j ' 

-1(a) 



It follows by (23) that there is a zero of L m m+J in the interval (0,C"i) and a zero in the interval 
(Cf-i i_i j Cj'i) f° r every i — 2, . . . , j . An analogous argument places zeros of L l ^l n+ - in the intervals 
( — Cm 1)0) an d ( — C m %t ~ C m -i i = 2, ...,m. By exhaustion, each of the above intervals 

contains one simple zero of L m " m +j- ^ 

We now study the distribution of the zeros of L m Q 2 as n — > 00. To that end, we will use the 
classical Heine-Mehler formula for Laguerre polynomials 

(24) (z/n) n- a =4 Z - a ' 2 J a {2^), n^oo, 

where J a (z) denotes the Bessel function of the first kind of order a (a > — 1) and the double arrow 
denotes uniform convergence in compact domains of the complex plane. 

The exceptional Laguerre polynomials admit a generalization of the classical Heine-Mehler for- 
mula, given by the following: 

Proposition 3.3 (Generalized Heine-Mehler formula). We have 



n —¥ 00. 



(25) L^(z/n)n- a ^(^ + 2 ^ J a (2^), 

A numerical representation of the convergence of the scaled exceptional Laguerre polynomials 
to the Bessel function is given in Figure [3] 



Proof. Multiply (21 ) by j a and replace z — > z/j. Taking the limit j — > 00 and using the classical 
Heine-Mehler formula (24) leads to 

J^Llt^ =3 (£a,m(0) -Urn-l(0))z- a/2 J a (2^)- 

covered by noting thi 

(0)-£a,m-l(0) 



The final expression ( 25 ) is recovered by noting that 

'a + m 
rn 
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as implied by (JTOJ) and ((22J). □ 

Note that for m = the classical Heine-Mehler formula is recovered as a particular case. 





Figure 1. Plot of n^L^^x/n) for m = 3, a = 5.5 and 20 < n < 100. The 
dashed red line corresponds to the limiting Bessel function x~ a ^ 2 ( a+J ^~ ") J a (2i/x) 
predicted by the generalized Heine-Mehler formula (25 1. 



Corollary 3.1. Let \zf \ be the sequence of zeros of the Bessel function J a {z) listed in 

I J i>l 

increasing order and let i denote the re 

Then we get the following asymptotic behaviour 



increasing order and let denote the regular zeros of L^^ n+ j(z) in the interval [0, oo). 



(26) 



lim jxj 



Proof. The above result follows from (251 and Hurwitz's theorem 



□ 



We have already seen that for fixed m the asymptotic behaviour of the regular zeros of the 
exceptional Laguerre polynomials coincides with that of the classical Laguerre. We now investigate 
in the same limit the behaviour of the m exceptional zeros of L I J^*l l+ j(z). 



i(fe) 

771,771+,?' 



in (—oo,0) converge to the m zeros of 



Proposition 3.4. As j — > oo the m zeros of L 
L^- l \-z). 

Proof. The classical Laguerre polynomials have the following outer ratio asymptotics 



(27) 



d a \z) 



Hence, by (21 ), we have for z £ [0, oo), 

-1(a) 



z ^ [0, oo), n — > oo. 



-(«) 



T l W TV 

m,m+j _ f. _ * j-1 f _ +. _ £ 

SCK,7Tl SO, 771— 1 _ 1 SCK,m SCK, 777— 1 S»— 1,777 - 



L 



(a) 



L 



Therefore, by Hurwitz's theorem the m exceptional zeros of ij^y converge to the zeros of 



l.m 



-(«-!) 



(-z). 



□ 
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Figure 2. Left: Exceptional zeros of the polynomials i^Ji+jC 2 ) f° r m = 6, 
a = 3.5 and 1 < j < 22. The squares denote the zeros of L„ 2) to which 
the zeros of L 1 ^ l ^ n+ -{z) converge for j — > 00. Right: Convergence is better seen by 
plotting the position of the negative zeros as a function of j. Solid lines show the 
limiting values. 



4. Type II Exceptional Laguerre polynomials 

4.1. Definition and identities. Let m > be an integer and a a real number. Let us introduce 
the polynomials 



(28) 



r) a ,m(z) ■= L ( m a) (z). 



For a fixed non-negative integer m and a real number a let us define the following first and 
second order operators 



(29) 
(30) 
(31) 
(32) 



+ (a - m)ri a+ i >m y, 

B a, m [y] ■= (v' -y)/va, m , 

C[y] := zy" + (a + 1 - z)y' , 

£a, m [y} : = £[y\ + 2z(iogVa+i, m )'(y - y') - my. 



The following factorizations follow from standard Laguerre identities: 

(33) C a = B^ m A^ m + a-m, 

(34) C l l m = ^J+i.rX+i,™ + a + 1 - m. 

For a real number a and given integers n > m > 0, we define the n th degree type II exceptional 
Laguerre polynomial by 



(35) 



L 



11(a) 



,11 rr( Q+1 )i 



j = n — m > 0. 



Expanding (29 1 and applying standard identities, the following dual representations of the type II 

l)Lt a - 2) Li a+1) 



polynomials hold: 

(36) Ljg^j = zL£*-»l£?> + ( 

(37) 



m — a 



= -zL^LW -(a + l+j)L^Lf\ 



Hence, up to a constant factor, the type II polynomials extend their classical counterparts, which 
are recovered for the particular case m = 0: 



(38) 



41 a) = -(l + « + n)LW. 
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The factorizations (33) and (34) yield the following intertwining relations between the standard 



Laguerre operator C a and the type II X-Laguerre operator r 
(40) B a+lrn C am = C a j r \B a+lm . 

The former relation provides the eigenvalue relation for the type II exceptional Laguerre polyno- 
mials: 



(41) 



-11(a) 
J m,m-\-j 



j = 0,1,2. 



The latter gives the following "lowering" relation between the type II exceptional Laguerre poly- 
nomials and their classical counterparts: 

(42) L^' = {l + n-a)L\- a -^Lf +1 \ j = n-m>0 

In order to find raising and lowering relations for the exceptional Laguerre polynomials, let us 
introduce the following first order linear differential operators 

Va+l, 
Va+l, 



(43) 
(44) 



B l lm[y] 



la+2, 



(y' - (logr?a+2,ro)'z/) , 

(zy' + (a + 1 - z)y) - z(\og-q a+2 , m )' 'y. 



In terms of these operators we have the following shape-invariant factorizations 



-11 _ All 



R n A 1 

a.m a.m ' 



F 

a. ■; 



B 

a,m a.m 



+ 1. 



(45) 
(46) 

From these factorizations the following lowering and raising relations for the exceptional polyno- 
mials easily follow: 

(47) A u |l ii(q) 

V ' J a.m m.n 



- T U( - a+1) n > m + 1 



(48) 



Bi 



jll(a+l) 



(n — m)L 



11(a) 
m,n+l ' 



n > m. 



The above equations can be conveniently re-written as 



(49) 



(50) 



L 



11(a) 



Va+2,m 



1]a-\-l,m 
Va+2.m 



L 



II(a+l) 



L m {a n +1) ) =(n-m + l) 



Va+2, 



Va+l,m / Va+2,m 



e Z r 11(a) 



4.2. Orthogonality. The type II exceptional Laguerre polynomials are formally orthogonal with 
respect to the weight 



(51) 



2 

Va+l.r. 



L 



(-a-1) 



The above weight is the solution W u = W of Pearson's equation ([8]) where 

p = z, q = (1 + a - z) - 2z(log r? a +i,m)' 



are extracted from (|32|). As a consequence, (|41|) and Green's formula imply 
(52) 



r 11(a) t- 11(a) 
m,ri2 ' m,rii 



fn 2 - m W L ii(q) W 11 = zVF 11 Wr 

where 712 > ni > m and where 

Wr[/, <?] = /'«?-/</ 

denotes the usual Wronskian operator. The following crucial result is established in 43 Ch. 6.73]. 
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Proposition 4.1. For a > m — 1 the polynomials r] a m (z) have no zeros in [0, oo). The number 
of negative real zeros is either or 1 according to whether m is even or odd, respectively. 

Thus, assuming a > m — 1 and restricting the interval of orthogonality to [0, oo), W^ m is a 



weight with finite moments of all orders, and the RHS of ( 52 ) vanishes, whch ensures genuine 
orthogonality in the I? sense. 

4.3. Zeros of the type II Laguerre polynomials. Henceforth, let us assume that a > m — 1, 
where m > is an integer. As above, we will call the real positive zeros of Lm'"'(z), n > m regular 



and the negative and complex zeros exceptional. From ( 36 1 we have 
r ii(a) / n x / ,fa + j + l\fm-a-i 

( 53 ) L m,m+j( ) = ( m + l )\ < j ){ m + l 

Hence, z = is never a zero of such a polynomial. 
Proposition 4.2. The zeros o/L^"'(z), n > ra, are simple. 



i = 0,l,2,. 



Proof. This follows by ( 36 ) ( 37 1 ( 42 ) and the fact that the zeros of the classical Laguerre polyno- 
mials are simple. 

Proposition 4.3. The polynomial L^n (2), n> m, has exactly n — m regular zeros. 



□ 



Proof. We prove the existence of at least j = n — m regular zeros by induction on j. The case 
j = is trivial. Suppose now that the proposition has been established for j > and a > m — 1. 
Since a + 1 > m — 1, the proposition is also true for Lm^ +1 '. Let Ci> • • • > Cfc> k > j,be the regular 
zeros of Lm[n +1 \z). By ( |50| ) and Rolle's theorem, I^^z) has at least one zero in each of the 
intervals (Ci, C2), ((2, Cs), ■ • • , (Cfc-i> Cfc)> Also by (50]), there is a zero in (0, Ci) and a zero in 00), 
for a total of at least A; + 1 zeros. 

We conclude by showing that j is also an upper bound for the number of regular zeros. The 
proof is again by induction on j = n — m. By (36), 

(54) = (m - 1 - a)^-"- 1 ) = (m - 1 - a)^ + i, m . 

by Proposition 4.1 the latter has no real, non- negative zeros. The lowering relation (49) shows 



that between two regular zeros of lS' at least one zero of n-i ( z ) nes - Hence, if we assume 
that the latter has at most j — 1 regular zeros, then the former has at most j regular zeros. □ 

Proposition 4.4. The type II polynomial L^, n> m, has either or 1 negative zeros, according 
to whether m is even or odd. 



Proof. Let 



m even 

1 m odd. 



Let e a . m ,n be the number of negative zeros of lll'n'- We wish to show that 

e Q ,m,n =e m , at > m - 1, n > m. 



Suppose that m is odd. By ( 53 ) 



sgn J L^>(0) = (-ir+ 1 . 



By (36) 



L 



Ufa) 



to — 1 — j — a 
m\j\ 



(—z) m+J + lower degree terms. 



As a consequence, 



lim sgnL^)(z) = -1. 
Hence, e a ,m,n > 1 if to is odd, and therefore e ai m in > e m . 



4.3 



L m°m,+j nas 3 + £a,m,n real zeros. Hence, by (|49|, i^m"+i-i has at least 

□ 



-II(a+l) 



( 54 1 and Proposition 



4.1 



By Proposition 

j — 1 + ( a ,m,n real zeros. Continuing inductively, £m,m~ has at least e Q!m ,„ real zeros. Hence, by 



£a,m,n < £m, as was to be shown. 



For the type II exceptional Laguerre polynomials, a Heine-Mehler type formula also holds: 
Proposition 4.5. As n — > oo, we have 

(55) - " ' " ' 



\ rn J 



Proof. Multiplying (37) by j 1 a , replacing z — > z/j, and applying the classical Heine-Mehler 
formula, we get 



The polynomials Lmm+j anc ^ ^™ " 1 ^ are uniformly continuous on compact subsets of C. Setting 
n = m + j, we have by uniform continuity on compact subsets 



4r a_1) (*A0 



m — a — 1 
m 



.— 1-a r n ( a ) 



as n — > oo. Equation ([9]) is needed to establish the first statement. 



□ 



Note that, as a consequence of (38 1, the above assertion reduces to the classical Heine-Mehler 
formula for m = 0. 

Corollary 4.1. Let < z\ < z 2 < z 3 < ■ ■ ■ denote the positive zeros of the Bessel function of 
the first kind J a (z) arranged in an increasing order and let < z n< \ < z n 2 < ■ ■ ■ < z nm _ n be the 
regular zeros of L^n also arranged in increasing order. Then, 



(56) 



lim 



^ 2 /4 



Proof. The above result follows from (551 and Hurwitz's theorem. 



□ 



Away from the interval of orthogonality, we can describe the asymptotic behaviour as follows: 
Proposition 4.6. As j — > oo we have 



- 11(a) 



;(*) 



L) >(z) 

on compact subsets of C/[0, oo). 

Proof. For the outer ratio asymptotics of the classical Laguerre polynomials, we have 



(-*) 



t/2 j 



L[ a+t \z) 



0(j*/ 2 ), i^oo 



uniformly on compact subsets of C/[0, oo). The desired conclusion now follows by (37). 



□ 



Proposition 4.7. As n — > oo the exceptional zeros of 1^ , n > m, converge to the zeros of 



.(-a-1) 



(z). 



Proof. The desired conclusion follows by the preceding Proposition and by Hurwitz's theorem. □ 
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Figure 3. Exceptional zeros of the polynomials L^^ + j(z) for m = 15, a = 14.01 
and 1 < j < 22. The squares denote the zeros of Lm a ^\z) to which the zeros of 
im,j +J (z) converge for j -> oo. 



5. Exceptional Jacobi polynomials 
5.1. Definitions and identities. Let m > be a fixed integer, and a,j3 real numbers. Let 

(57) T a j[y] = (1 - z 2 )y" + (/?-« + ( a + / 3 + 2)z)y', 

denote the Jacobi differential operator. The classical Jacobi polynomial of degree n can be defined 
as the polynomial solution y — P r [ a '^\z) of the second order linear differential equation 

m - r-.,_ ./ a , _(« + !), 

Next, define 
(59) 

(60) 
(61) 



r Q ,^[y] = -n(l + a + /3 + n)y, y{l) 

u 

T a ,p,m[y] = T at p[y] +(a-/3-m+ l)my 

- (logP^- 1 ^- 1 ))' (/?(1 - z)y + (1 - z 2 )y') , 
A a , , m [y] - (1 - ^)Pir Q '« y' + (m - aJP^" 01 - 1 ^- 1 ) y, 

r> r -I (1 + + (1 + 



P 

The following operator factorizations can be verified by the application of elementary identities. 

(62) T a ^ <m = -A 0+ i )J 8_i im .B a+ i j( 8_i i7n - (m — a — l)(m + 

(63) T a ^ = B a ^ >m A a> p tTn - (m - a)(m + £ + 1). 
For n > m, we define the degree n exceptional Jacobi polynomial to be 



(64) 
(65) 



p(«,/3) = 



a + 1 + j 

J-l) 
a 



l,/3-l,m 



P 



j = n — m > 0, 



^ Q(l + a + /J + i)(z - l)Pi— Wpp** 



+{a + 1 - m) P^ 



(-o-2,/3)p(a+l,/3-l) 



13 



p("./3) _ p(a,/3) 
r 0,n ~ r n 



The exceptional polynomials and operator extend their classical counterparts 

(66) T a,p,o[y] = T <*Ay}> 

(67) 

By construction, these polynomials satisfy several identities, which we enumerate below. The 
factorizations ( 62 ) ( 63 ) give the intertwining relations 

(68) Aa+xfi-x^mTa-it-ifi-i — T a ^^ m A a+ i^p^i^ m , 

(69) T a+ i^_iB a+ i^_i^ m = B a+ i 

,{3 — 1,771 

From the above relations we can derive the eigenvalue equation for the X m -Jacobi polynomials 



(70) 



P£f \ = ~(n - m)(l + a + p + n - m)P£f , n>m. 
The factorization ( 63 ) implies the following identity 

(71) M)"> + 1 + j)(0P^ +j + + l)Pjn,£.s) = 

(a + 1 - m + j)((3 + m + fi P <r-«-i,P-D pf+W-Q . 



It will be useful to express P in a way that is symmetric in the dimension j and the codimension 
m. Namely, 



(72) 
(73) 



(-ir(a+i)P£^ +1) 

= (a - m)P^ a - 1 ^P^ + (z- l)Pt a ^P^ Y 
= (a + i)Pf ^M"^ - (z- l)P,<-»J&-.0'. 



The first equation is just a restatement of the definition (64), while the second identity follows 
from the classical relation 

(74) (z - l)pf = aPf'« - (a + j)pf " W) . 

At the endpoints of the interval of orthogonality we have the following classical identities 

(75) P^)(-l) = (-!)" + 

n! 

(76) pM {1) 

n! 

which in the case of exceptional Jacobi polynomials yield the following generalizations 



(77) 
(78) 
(79) 



J m,n V- 1 -/ 



a + n — m 
n 



(a + 1 - m) m+ j 



m! j\ 



n > m, 



j = n — hi, 



A(a,« / lV - (^ + j + m)(l + a-m + j)()3 + l) m -i(/3) 3 - 

m ' m+jl j 1 J (l + a + j) m!j! 

Define the lst-order operators 



(80) 



(81) 



-4 



-2,/3) 



a,/3,m 



B, 



p 

1 , 



± 7 



(-a-1,/3-1) r 



a,P,mM = (l-Z )— 

' f f 



/3) 



y'- (logP4-«-^- 1 ))' + 



a + 1 )8 + l 



1 - z 



1 + z 
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The following "shape-invariant" factorizations relate exceptional operators of the same codimcnsion 
at different values of the parameters a, P 

(82) 

(83) A a ^^ m B a ^^ m = T a+ i^ + ij n + 2 + a + j3. 
The corresponding intertwining relations, namely, 

(84) T a+ i t p + i im A 

(85) Bo,ftmT a +l^+l : m = Ta,/3,m-Sa,0,m, 

give rise to the lowering and raising relations for the exceptional Jacobi polynomials 

/ p{a,p) \' 1 / p (- a -l,p-l)\ 2 p( a +lJ3+l) 

( 86 ) pC-a'-a.g) - g(n-m + a + /3 + l) "' ( _ a _ 2 ^ (-^-i) ' ri>m + l, 

(87) (1 - ,)-(! + «)-/> (1 - + z)^ +1 4^ 



p(-a-l,/?-l) 



/ pi-a-2,0) Y p(a,/3) 
o/ 1 , -i \ / - r ™ 1 m.n+1 ^ 

-2(n-m + l) ( _ a _ 2</?) , u>m. 

\ ' m / J m 

As usual, we denote by f'(z) the derivative of / with respect to the z variable. 

5.2. Orthogonality. The exceptional Jacobi polynomials are formally orthogonal with respect to 
W a fi, m (z)dz, — 1 < z < 1, where 

(88) W a , fi , m {z) = r (1 ~ z)a(1 + ^ 2 . 

In order to have orthogonality in the L 2 sense, additional conditions need to be imposed on the 
parameters a,/3, and m. The condition a,P > —1 is necessary for the measure (88) to have finite 
moments of all orders. Another requirement is that the denominator ?i a 1,13 ^does not vanish 
for z 6 (—1, 1), which imposes extra conditions on a, /3, and m. 

An analysis of the zeros of classical Jacobi polynomials can be found in Szego's book |43l Chapter 
6.72] . First, let us recall that P^'^ (z) has a zero of multiplicity k at z = 1 if a = —k, k = 1, . . . , m, 
and a zero of multiplicity j at z — — 1 if p = —j, j = 1, . . . , m. 

We also mention the degenerate cases where 

degP£* ,/3) = fc when n + a + (3 + l = -k, k = 0, 1, . . . , n - 1. 

For such parameter values the nth Jacobi polynomial has degree < n. In these degenerate cases 
(where a + (i is a negative integer) , we have 

(89) ( a + m )p^\z)=( a + n )p^\z), a + P = -l-m-n. 
\ m J \ n J 

Since f3 > — 1, the denominator has a zero at z = — 1 if and only if p = 0. However, the latter 
condition gives a weight with an overall factor of (1 + z) , which would violate the assumption 
that it has finite moments of all orders. Therefore we must impose /? ^ 0. The condition that 
Pm a ^ 1 ' l3 ^ 1 \z) 7^ for z G (—1, 1) is satisfied in exactly two cases 

(A) Both P and a + 1 - m e (-1, 0). 

(B) Both P and a + 1 - m G (0, oo). 

For (A) we have m — 2 < a < m — 1, while for (B) we have a > m — 1. Therefore, in both cases 
Pir 1 " Q,/3 " 1) (l) 7^ 0. According to identity Q, we also require 

(90) a+l-m- P {0, 1, . . . , m - 1}. 
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If this condition is violated, then degP™ a 1,13 X \z) < m and, therefore, the codimension (see 



below for discussion) is a + 1 — m — (3, rather than m. We therefore append condition (90) to the 
assumptions in (A) and (B), to complete the following 

Proposition 5.1. Suppose that m > 1. The measure W a> p tm dz, z 6 (—1,1) is positive definite 
with finite moments of all orders if and only if a, /?, m satisfy one of the following conditions 

(A) p,a+l-me (-1,0). 

(B) p,a+l-me (0,oo). 

In order to ensure that degPm a ^ — m it is also necessary to require that a + 1 — m — fj ^ 
{0,l,...,m-l}. 

5.3. Exceptional Flag. Let us define the following codimension m polynomial flag {U a) p tm J\JLy 
where 

U a ,p, m , 3 = {f{z) G P m+j -i(z) : P^ a -^-V\(l + z)f'+/3f}- 



At the level of flags, the factorizations (62) (63) correspond to the linear isomorphisms 



a+l,P-l,m '■ "Pj-l E/o,ftm,ji Ba+l,^-l,m ■ U a ,/3,mJ ~> Pj-1> j — 1)2, 

Thus, the exceptional polynomials P^m+j gi ye a basis of the flag {?7o,^,mj}|^i 



Proposition 5.2. Suppose that a, /3, m satisfy either condition (A) or condition (B). Then Pmm+j( z ) 
is the unique polynomial in U a p m j, orthogonal to !7 a ,8,mj-i with respect to W a « m dz, z G 



(—1, 1) that satisfies the normalizing condition (77) 



Once we have analyzed the underlying factorizations that give rise to A m -Jacobi polynomials, 
and the conditions on the parameters that ensure their L 2 -orthogonality, we can now turn our 
attention to describing some properties of their zeros. 

5.4. Zeros of exceptional Jacobi polynomials. Let us refer to the real zeros of Pmm+j( z ) m 
z G (— 1, 1) as the regular zeros. All other zeros, whether in (— oo, —1) U (l,oo), or complex, will 
be said to be exceptional zeros. 

Proposition 5.3. Suppose that a,f3,m obey either condition (A) or condition (B). Then the 
regular zeros of P^'m+j ( z ) are simple. 



Proof. This follows from (64) (71) and the simplicity of the zeros of classical Jacobi polynomials. 

□ 

Proposition 5.4. Suppose that a,/3,m obey either condition (A) or condition (B). Then Pm 'm+j 
has exactly j regular zeros and m exceptional zeros. 

Proof. We begin by showing that P^'m+j has at least j regular zeros by induction on j. The case 
j = is trivial. Suppose now that the proposition has been established for j > 0. Note that 



a + 1,(3 + l,m always belong to class B by hypothesis. We observe also 43 Chapter 6.72] that 
Pi - "' 1 '^ X \ z ) and Pm a 2 ' l3 \z) have no zeros in z G [—1, 1]. Let Cl> •• • iCti i>j, be the regular 



zeros of Pm,m+j +1 ^ ■ By (87) and Rolle's theorem, Pm,m+j+i( z ) has at least one zero in each of the 
intervals (Ci, C2), (C2 , C3); • • • j (Ci-i> Ci)- There will also be a zero in (—1, Ci) and a zero in (£j, 1), 
for a total of at least i + 1 zeros. 

We conclude by showing that Pj^ 'm+j has at most j regular zeros. The proof is again by 
induction on j. Observe that 

a + 1 — m 
a + 1 + j 



p(<x,P) — f_1 \rn I 



p(-a- 



2. 



and the latter has no zeros in [—1,1]. Relation (86) shows that between two regular zeros of 



Pjnm+j there is at least one zero of Pm^n+j-i- H ence > if we assume that the latter has at most 
j — 1 regular zeros, then the former has at most j regular zeros. □ 
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5.5. Asymptotic behaviour of the zeros. Our next goal is to derive a representation for the 
X m -Jacobi polynomials that is amenable to asymptotic analysis. 

Proposition 5.5. The following identity holds: 



(9i) i-irp^ilj 



a - (3 - m + 1 p (- a ,p) 



J 



-P 



1 + a+j ' m ~ x \a + P + If 1 a + P + 2j i_1 

1 + q-m pf-a-qjn , 3 p (- a - h0 -i)\ p ( a ,p) 
1 + a + j m ^ 1 + a+j m 1 ' 



+ 



Proof. We begin with (73 1 and apply (74) as well as the following classical identities: 

p(l+oi,/3-l) _ p {a+l,p) + p {<x,0) 

(a + /3 + 2j)(z l)P^t m = 2jPj°" P) - 2(a + j)p£f 



□ 



Proposition 5.6. Suppose that a, ft > — 1. ^4s j — > oo we ftawe the following asymptotic behaviour 
for z in compact sets ofC\[— 1, 1] 



(92) 



6(o./3) 



DC'. 



Proof. We make use of the following well known ratio asymptotics formula for classical Jacobi 
polynomials: 



(93) 



»(z) 



1 



P 



**[-M]. 



The conclusion now follows directly from (91). 



□ 



As a straightforward consequence, the following corollary describes the asymptotic behaviour 
of the zeros of exceptional Jacobi polynomials. 

Corollary 5.1. The j regular zeros o/P^m+j approach the zeros of the classical Jacobi polynomial 
Pj 01 as j oo, while the m exceptional zeros of P^'J^ + j approach the zeros of Pm a 1 \ 

The Heine-Mehler formula for the classical Jacobi polynomials states 

(94) n- a P 7 [ a ^{cos(z/n)) =t (z/2y a J a (z). 

The X m -Jacobi polynomials satisfy a generalized Heine-Mehler formula, given by the following 
proposition. 



Proposition 5.7. When n — > oo, we get 

(95) n- a P£f (oos(Vn)) =4 (-l) m (z/2)- a f 



m — 1 — a 
m 



J a (z), n> m. 



Proof. Taking the limit j — > oo in (91 ) and using the classical Heine-Mehler formula (94) leads to 



the desired result, keeping in mind also that 

Pi"^- 1 )(l)=( 



m — 1 — a 
m 



a 
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6. Summary and Open problems 

We have provided suitable representations of exceptional polynomials in terms of their classical 
counterparts by exploiting the isospectral Darboux transformations that connect them. These rep- 
resentations allow to derive Heine-Mehler type formulas for the exceptional Jacobi and Laguerre 
polynomials, which describe the asymptotic behaviour of their regular zeros (those lying in the 
interval of orthogonality). The behaviour of the regular zeros of exceptional polynomials follows 
the same Bessel asymptotics as the zeros of their classical counterparts. We have also proved 
interlacing between the zeros of exceptional and classical polynomials, while the zeros of consec- 
utive exceptional polynomials also interlace according to their Sturm-Liouville character. As for 
the exceptional zeros (those lying outside the interval of orthogonality) we have established their 
number and location and we have proved that for fixed codimension m and large degree n they 
approach the zeros of a classical polynomial. We have performed a careful analysis of the admis- 
sible ranges of the parameters that ensure a well defined Sturm-Liouville problem. We have also 
given raising and lowering relations for the exceptional polynomials. These relations correspond 
to a shape-invariant factorization, i.e. a Darboux transformation that falls within the same class 



of operators, with a shift in the parameters (see for instance ( 43 1— ( 46 )) , and they imply that the 
associated potentials in quantum mechanics will be exactly solvable and shape invariant. 

It was recently noticed that more families of exceptional orthogonal polynomials can be con- 
structed through multi-step Darboux or Darboux-Crum transformations [18], an idea that has 



been further developed in 23 35 39 . In this work we have analyzed the zeros of exceptional 
orthogonal polynomials that can be obtained from the classical ones by a 1-step Darboux trans- 
formation. These polynomials can be written as a first order linear differential operator acting on 
their classical counterparts and the exceptional weight is a classical weight divided by the square 
of a classical polynomial with zeros outside the interval of orthogonality. An open problem is to 
extend this analysis to multi-step exceptional families, where exceptional polynomials are obtained 
by the action of an m-th order differential operator. 

We believe that all exceptional orthogonal polynomials can be obtained from a classical system 



by a multi-step Darboux transformation, 20 , and the exceptional weight for these systems will 
have in its denominator the Wronskian of all the factorizing functions, which are essentially classical 
orthogonal polynomials. The characterization of all such Wronskians whose zeros lie outside the 
interval of orthogonality becomes then a crucial question. 

It is trivial to know the location of the zeros of a classical polynomial, and therefore to constrain 
its parameters so that they fall outside the interval of orthogonality for the exceptional weight to 
be regular. The question becomes much more involved when dealing with a Wronskian of classical 
polynomials as it happens in the multi-step case. However, this question must be addressed in 
order to select those multistep weights that are non-singular. 

The position of the zeros for Wronskians of consecutive Hermite polynomials have been inves- 
tigated numerically by Clarkson [4] since these functions appear as rational solutions to nonlinear 
differential equations of Painleve type, j5j. A further numerical analysis together with some con- 
jectures in a more general case have been recently put forward by Felder et al. 10 , in connection 
with the theory of monodromy free potentials. We stress that a Wronskian of classical polynomials 
might have no real zeros even if the polynomials themselves do. The Adler-Krein theorem 



provides a useful criterion to identify these cases, and it is actually a much more general result 
for eigenfunctions of a Schrodinger operator, not just polynomials. A generalization of this re- 
sult is being carried out by Grandati, who is extending the analysis to factorizing functions of 



isospectral Darboux transformations 23 , as opposed to the Adler-Krein case which refers only to 
state-deleting Darboux transformations, for which the factorizing functions are true L 2 eigenfunc- 
tions. The most general problem of Wronskians that involve factorizing functions of mixed type 
remains unsolved. 
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